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Abstract 

We use retrodictive quantum theory to describe cavity field measurements by 
successive atomic detections in the micromaser. We calculate the state of the 
micromaser cavity field prior to detection of sequences of atoms in either the 
excited or ground state, for atoms which are initially prepared in the excited 
state. This provides the POM elements which describe such sequences of 
measurements. 
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I. INTRODUCTION 



The advent of quantum information theory, coupled with an increased technological 
ability to produce and study exotic entangled states of the optical field has stimulated a 
re-examination of the fundamental principles of quantum mechanics. One example of this is 
the recent resurgence of interest in retrodictive quantum theory, in which the likelihood of 
past events is inferred from measurements of the present. This subject was, until relatively 
recently, a simple problem associated with time-reversal in quantum theory UJ. It has now 
evolved into a powerful technique for calculating retrodictive conditional probabilities in 
quantum theory, particularly when the knowledge of the output state produced by a state 
preparation device is not complete. The crucial step in this was the use of Bayes' theorem 
|| coupled with normal predictive quantum theory. This allows a retrodictive state to be 
defined such that retrodictive and predictive conditional probabilities are consistent ||. The 
theory can be applied to both closed systems, in which the full evolution can be followed, 
and to open ones, in which the system of interest interacts with an unmeasured environment 

The predictive state of a system at any time during the evolution is the prepared state 
evolved forwards in time. The 'collapse of the wavefunction' occurs at the measurement 
device when this evolved state is projected onto the measurement basis. The retrodictive 
state at any time between preparation and measurement is the measured state evolved 
backwards in time. At the preparation time this state collapses onto any one of a set of 
possible initially prepared states. In open systems, where the system must be described 
by a density operator, we can define a retrodictive density operator which can be evolved 
backwards in time according to a retrodictive master equation |J. 

The micromaser has been useful tool for quantum optical and cavity QED experiments 
for many years ||. In the basic device highly-excited Rydberg atoms are passed through a 
high-Q microwave cavity, one of whose longitudinal modes has a frequency which is close to 
one of the atomic transition frequencies. The atoms interact with the cavity field, and the 
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whole device behaves to a very good approximation like a two-level atom interacting with 
a single-mode quantum optical field. The atoms can be prepared and measured in what 
amounts to either the excited or ground state. In a further refinement resonant microwave 
pulses can be applied to the atoms prior to entering and after leaving the cavity. This allows 
the atoms to be prepared or measured in any desired superposition of the excited or ground 
states. The micromaser operates in what has been termed the strong-coupling regime of 
cavity quantum electrodynamics, in which a single photon in the cavity is a strong enough 
field to cause multiple excitations and de-excitations of the atom. 

Usually there is no direct way to measure the cavity field in the micromaser. The only 
way to tell whether or not the cavity field has been excited by an atom is to send a further 
atom into the cavity in a known state and to measure the atomic state after it has passed 
through the cavity. It would therefore be useful to know precisely what measurement an 
atomic detection in a particular state (or a sequence of detections in various states) makes 
on the cavity field. This paper is designed to solve this problem for the case that all of the 
atoms enter the cavity in the excited state. If an atom or sequence of atoms in a known 
state are passed through the cavity and detected as a sequence of atoms at the output, what 
does this particular sequence of detections say about the cavity field prior to the passage of 
the atoms? 

This paper is organised as follows. In Section (II) we give a brief overview of retrodictive 
quantum theory. Section (III) describes how this may be applied to the two-level atom 
interacting with a single field mode, with particular emphasis on retrodicting micromaser 
field states form atomic state measurements. Section (IV) discusses the results of particular 
measurements and in Section (V) contains a summary and discussion of the results. 

II. QUANTUM RETRODICTION 

In this section we provide brief details of retrodictive quantum theory. A fuller account 
can be found in references Elsewhere in this volume we provide a more complete 
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discussion of the quantum theory of preparation and measurement ||. The analysis which 
follows, based on preparation operators and positive operator measures (POMs) (T(J is nec- 
essarily incomplete, but sufficient for our purposes. 

Suppose that we have a preparation device which produces output states pf re with prior 
probabilities P{i)- pf red is the usual density operator of predictive quantum mechanics. 
This state can evolve and interact with other systems until it is measured by a measuring 
device. A measurement in which none of the measurement results are lost or discarded can 
be described by a measurement POM [ 10| . This is a set of non-negative definite, Hermitian 



elements flj which sum to the unit operator, each element corresponding to a particular 
measurement outcome. In general there is no requirement that there be the same number 
of POM elements as there are states which span the system space, but for von Neumann 
measurement this is so, and the POM elements are simply the projectors of the particular 
chosen states which span the space. Suppose that preparation takes place at time tp and 
measurement at a later time t m . Within this framework the predictive probability that the 
measurement outcome flj is obtained given that the state pf red was prepared is 

P(j|.) = Tr(pf red (t m )fl,), (1) 

where 

pf red (tm) = U(r)pr%)UHr) (2) 

is the initial density operator evolved forwards in time to the collapse time and 

U(t) = exp (3) 

is the evolution operator, which operates for the length of time between preparation and 
measurement, r = t m — tr>. 

Suppose that instead of calculating the predictive probability P(j\i) we wish to calculate 
the retrodictive conditional probability P(i\j) that the state pf ve was prepared, given 
a particular measurement result IT,. It is possible to do this by calculating all possible 



predictive conditional probabilities for the system, and then using Bayes' theorem. A simpler 
and more natural approach is to use retrodictive quantum theory, so that the required 
probability can be written 



Tr 



A*pf tr (*p) 



P( ' |j) Tr [Apf tr (i p )] ' (4) 
Here the operator Aj is the preparation device operator, and 

A = EA i = E^)pF red (5) 

i i 

is the a priori density operator, the sum of each possible preparation density operator 
weighted by its prior probability of production. A is the best description of the state we 
can give without knowing the outcome of the preparation or measurement. The retrodic- 
tive density operator at the preparation time is simply the normalised measurement POM 
element evolved back from the measurement time to the preparation time, 

P? *%) = U\r)pf(t m )U(r), (6) 

with 

= (7) 

The above formulae for the conditional probabilities (eqs. ([!]) and @) apply equally well 
for open systems, where the system of interest interacts with an unmeasured environment 
with many degrees of freedom. If this environment causes information to be lost about the 
system, and the Born-Markov approximation holds, then the the evolved density operators 
are the solutions of master equations [|1TJ. In eq.(|I]) the density operator required for such 



an open system, which for a closed system would be given by eq. (0), is the solution of 
the usual master equation forwards in time from the preparation time to the measurement 
time. However, in eq. (f|) the solution required instead of eq. (|6|) is that of the retrodictive 
master equation, giving the evolution backwards in time from the measurement time to the 
preparation time. We have recently derived this equation from the general principle that 



conditional probabilities should be independent of the time of collapse of the wavefunction 
||. This principle, which applies equally well in open or closed systems, leads to an unusual 
interpretation of measurements in a quantum system. The measurement POM elements, 
when evolved back in time, still satisfy the criterion that they form a valid POM, namely 
that they sum to the unit operator. This means that the evolved POM still represents a 
valid measurement of the quantum system, and we may regard any amount of the evolution 
as still forming part of the measurement. The same may be said of the prepared density 
operator evolved forwards in time. 

In the system considered in the present paper, however, the Born-Markov approximation 
is not made, and so we must consider the full evolution of the coupled atom-field system. 



III. PREDICTION AND RETRODICTION FOR THE COUPLED ATOM-FIELD 

SYSTEM 

Here we apply the retrodictive formalism to a coupled system consisting of a two-level 
atom and a single cavity mode of an electromagnetic field. We have already applied retro- 
dictive quantum mechanics to this system. When the field is excited by a coherent state, 
predictively there are collapses and revivals of the Rabi oscillations in the atomic excitation 
probability. If the atom is measured to be in a particular state there also exist earlier revivals 
or 'previvals' in the atomic excitation at times prior to the measurement [[7|. Here we will 
concentrate specifically on the micromaser, in which atoms pass through a microwave cavity 
one at a time. The interaction between an atom with upper level |e) and lower level \g), 



and an electromagnetic field is governed by the Jaynes-Cummings Hamiltonian ||11|| . In the 
interaction picture this is 

H = ~^^3 ~ "^A (o + a — a^aJ^j , (8) 

where A is the detuning between the atomic frequency and the cavity mode, = |e)(e| — 
\g)(g\ is the atomic inversion operator, a + = \e)(g\ and <r_ = |<?)(e| are the atomic raising 
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and lowering operators, a) and a are the creation and annihilation operators for the single 
mode field, and A is the coupling constant. The rotating wave approximation, which has 
been made in deriving this Hamiltonian, ensures that whenever a photon is lost from the 
field the atomic state must change from \g) to |e). In the standard predictive picture of 
quantum mechanics a coupled atom-field system evolves forwards in time according to this 
Hamiltonian from a preparation time tp to a measurement time t m . After this has happened 
the coupled density operator for the whole system is 

plT d (t m ) = U(r)pf e %) ® pr%)U\r), (9) 

where p^ ed (tm) is the coupled density operator for the atom-field system at the measure- 
ment time, p a re ^(£p) an d Pj^^itp) are the uncoupled atom and field density operators at 
the preparation time. U(r) is the evolution operator for the coupled system, given by eq.([jD 
with Hamiltonian given by eq. @. In the micromaser we generally wish to infer the state 
of the cavity field from a measurement of the atom. This is done by conditioning on the 
atomic measurement outcome and tracing over the atomic states, 



pf ed (t m ) oc Tr a \u & U(r)pf e %) ® /f ^p)^) 



(10) 



where n a is the POM element corresponding to the outcome of the measurement performed 
on the atom. If no measurement is performed on the atom then fl a = i a . Similarly we may 
not have any information about the initial state of the field, in which case p^ re< ^(tp) oc If. 

The retrodictive picture differs from above in that the state of the system before the 
measurement is assigned on the basis of the measurement outcome. Thus if the measurement 
POM elements for the atom and the field are fl a (tm) and f[f(i m ), the coupled initial density 
operator corresponding to equation @ is 

Paf(*p) = ^)Pa etr (*m) ® pf* (t m )U (r) 

oc U\r)fl & (t m ) ® n f (t m )C/(T). (11) 

In the micromaser the atom will in general have been prepared in a particular initial state, 



and this state conditions the coupled density operator similarly to equation (fLOf), to give the 
retrodictive density operator for the field 

pp tr (tp) oc Tr a [pf ed (t p )^( r )n a ( tm ) g, fl f (t m ) f >( r )] , (12) 

where the constant of proportionality is determined by normalisation. If there is no informa- 
tion at all about the preparation of the initial states then Pa ied (£p) becomes proportional 
to the unit operator for the atomic space. Equation fll2|) is enough to provide all of the 
information needed for this paper. 

IV. RETRODICTION OF MICROMASER FIELD STATES 

In this section we apply the formalism derived in the previous section to the micromaser. 
Two-level atoms in their excited state are passed through an optical cavity and a measure- 
ment of their state is performed. We assume that the atomic frequency is resonant with 
the micromaser field mode. This is a typical set-up of a micromaser experiment [P,l2|. The 



final field state is unmeasured and we wish to know what the initial field state was before 
the atoms passed through. The formalism can equally well be applied for any other initial 
atomic state, but we will leave this for a further publication. In order to find the initial field 
state from a sequence of atomic measurements we first consider the effect of a single atom 
passing through the cavity. Equation fli~2|) provides the conditioned retrodictive field state 
given knowledge of the initial and final atomic states. The preparation time tp is effectively 
the time at which the atom enters the cavity, and as the atomic state does not change after 
the atom leaves the cavity, the exit time is the measurement time t m . 

Before their passage through the cavity the atoms are initially in their excited state, so 
Pa(^p) = |e) (e| . The trace over the atomic states is trivial, and the initial field state is 

pf tr (tp) ex (e\U\r)U a (t m ) ® U f (t m )U(r)\e). (13) 

We assume that the cavity POM element at t m is given by 
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n f (t m ) cc p f (t m ) = Y,Pn{tm)\n)(n\, (14) 

n 

where P n {tm) is the probability that the field contains n photons at (t m ). If the atom is 
detected in the excited state then the POM element (and retrodictive density operator) is 
simply the excited state projector. The unitary evolution operator is given by equation ([3|). 
A little algebra gives the result 

pf tr (t p ) oc^P n (t m )cos 2 [fi(n + l)r/2]|n)(n|, (15) 

n 

where f2(n) = 2\^/n. Thus the initial state photon number probabilities are equal to the 
final state probabilities multiplied by a cosine function whose argument is the product of 
the Rabi frequency and the interaction time. If the atom is detected in the ground state a 
similar calculation gives 

pp tr (tp) oc^P n+1 (t m )sin 2 [fi(n + l)r/2]|n)(n|. (16) 

n 

In this case the probability that the initial state contains n photons depends upon the final 
state probability for n + 1 photons and a sine function which contains the n + 1-photon Rabi 
frequency. The occurrence of P n +i(t m ) in the above formula is easy to explain, as when the 
atom exits in the ground state there is one more photon in the final field than the initial. 

Successive atoms entering the cavity can be treated one at a time. If two atoms enter 
the cavity in sequence then the initial retrodictive state into which the cavity is projected 
by the detection of the second atom becomes the final state in an identical calculation to 
the above for the first atom entering the cavity. Thus every time an atom enters the cavity 
and leaves in its excited state the photon number probabilities undergo the transformation 

P n (t v ) cc P n {t m ) cos 2 [n(n + l)r/2], (17) 

and when an atom leaves in the ground state the transition is 

P n (*p) ex P n+ i(t m ) sm 2 [n(n + l)r/2]. (18) 
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The exact retrodictive density operator can be found by normalisation at the end of the 
calculation by simply requiring all of the probabilities at tp for the first atom to sum to 
unity. 

As a simple example consider two atoms passing through the cavity. There are four 
possible pairs of detection events which are summarised in the table below. 



detection events 


initial photon number relative probability 


e) then 


e) 


J R„cos 4 [fi(n + l)r/2] 


\g) then \g) 


P n+2 sm 2 [Q(n + 2)r/2] sm 2 [Q(n + l)r/2] 


e) then \g) 


P n+1 cos 2 [Q(n + l)r/2] sin 2 [fi(n + l)r/2] 


\g) then 


e) 


P n+1 cos 2 [Q(n + 2)r/2] sin 2 [Q(n + l)r/2] 



Table 1. Relative probability that the cavity contains n photons prior to the passage and detection 
of two atoms. P n is the final state probability for n photons after both atoms have been detected. 



It is clear from the last two lines of the table that the order of detection matters. Every 
time an atom is detected in its ground state the number of photons in the cavity increases, 
and the Rabi frequencies increases correspondingly. Whether this occurs before or after the 
detection of an atom in the excited state will affect the Rabi factor imposed on the initial 
state probability. 



Measurement outcomes are characterised by their corresponding POM elements |nj • The 
atom detection experiments make a measurement of the cavity field, so there must be a POM 
element corresponding to this measurement. Different numbers and orderings of atomic 



detections in either the excited or ground states correspond to different POM elements [10 
For s atomic detections in either the excited or ground state there are a possible 2 s POM 
elements fti s . These POM elements can be used to calculate the probability for any given 
sequence of atom measurements given a known initial field state density operator pf(tp). 
The POM elements may be found by using an unmeasured state as the final cavity state 
that, in the usual predictive formalism, would have the unit operator for the cavity state 
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space as its associated POM element. This implies no knowledge of the final cavity state. 
It is convenient to set an upper limit N on the allowed photon number distribution. This 
allows us to determine a well-behaved retrodictive field state. A derivation of the field POM 
elements without this cut-off is given in the appendix. Thus in place of the unit operator 
we can use the operator 

N 

if=EI»X»l ( 19 ) 

n=0 

corresponding to a final state lf/(iV + 1). This is the unit operator acting on the iV + 
1-dimensional subspace spanned by the photon number states up to and including \N). 
Backward time evolution including the passage through the cavity of each of the s atoms 
yields the POM elements ITi. s in the form 

N 

ni.. s = £C nil .»(n| (20) 

n=0 

where the coefficients C n ^„ s will be proportional to the corresponding coefficients in the 
expression for pp^ r (tp). The POM elements must be normalised to sum to the unit operator 
or, in this case, the unit operator If that acts on the restricted photon subspace. For the two- 
atom case the coefficients C n ^.. s for the four POM elements are just the trigonometric factors 
given in Table 1. For example C n ^ ee for the detection of |e) then |e) is cos 4 [f2(n + l)r/2]. 

Figures 1-4 show photon number probabilities (and thus give the POM elements) corre- 
sponding to particular sequences of atomic detections. In each figure the interaction strength 
and time have been set so that At = ir. 

Figure 1 shows the relative photon number probabilities for a sequence of atoms detected 
in the excited state. If a single atom passes through the cavity, and is detected the initial 
photon number distribution of the cavity has the cosinusoidal form shown in figure 1(a). 
The distribution has maxima at photon numbers which are one less than perfect squares, 
corresponding to Rabi phase shifts of integral multiples of 7r. If the cavity contains one of 
these numbers of photons then an atom in the excited state must exit in the excited state. 
This is what is known as the trapping state condition, and it has been used to prepare 
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photon number states of the electromagnetic field ||12|| . If a whole series of atoms pass 
through and are detected in the excited state the troughs in the distribution become wider 
as the passage of one more atom means that each photon number probability gets multiplied 
by a cosinusoidal factor smaller than unity for all photon numbers except those one less than 
perfect squares. This is shown for five atoms in figure 1(b). 

Detecting atoms in the ground state gives rise to a different cavity POM element. In 
this case the minima of the sinusoidal distribution fall at photon numbers which are one less 
than perfect squares. One of the photons in the cavity when the atom is measured had to 
come from atomic de-excitation so there had to be one fewer photon in the cavity prior to 
the passage of the atom. Thus the whole distribution at t m is shifted by one photon at t p . 
The distribution for one atom detected in the ground state is shown in figure 2(a). 

Detecting atoms in the ground state thus has two effects: (1) it shifts the photon number 
probability distribution downwards by one photon for each atom detected and, (2) for a 
system which satisfies the trapping state condition it introduces zeros into the distribution . 
The zeros occur because the Rabi phase shifts for excitation numbers one greater than these 
photon numbers are tt, and so atoms prepared in the excited state can not exit in the ground 
state. These two effects combine to allow the measurement of photon number distributions 
with large gaps in them. For example, if three atoms are detected in the ground state then 
the cavity field cannot have contained fewer than four photons, as shown by figure 2(b). 
This POM element also has a photon number gap from 6-8 photons, and there are also 
gaps at higher photon number. These gaps extend their range downwards by one photon 
with each extra ground state detection. This is shown in figure 2(c), which is the photon 
number distribution prior to six atoms detected in the ground state. The initial cavity 
photon distribution must have contained at least nine photons. This also gives information 
about the number of photons in the cavity after the detection. There must be at least 15. 
Six of these come from the passing atoms and at least nine were in the cavity initially. 

If we have a little extra prior information about the cavity photon number distribution 
this measurement can allow us to measure and prepare specific photon number states. For 
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example, say two atoms are sent through the cavity and detected in the ground state. If we 
know that there are not more than five photons in the cavity after the atoms have passed 
through (not more than three prior to the atomic passage) then these two atoms are sufficient 
to give an initial cavity distribution of one photon for certain (figure 3). This measures the 
one photon number state. Furthermore, the fact that the cavity initially contained one 
photon means that the de-excitation of the two atoms has now created a three photon state. 
Other possible measurements and preparations are possible, with the restriction on photon 
number a little higher. 

In general a set of atoms passing through the cavity will not all be detected in the excited 
or ground state. Detection of any atoms in the ground state will introduce zeros into the 
photon number distribution, and these will move towards lower photon numbers with each 
successive atom detected in the ground state. Otherwise the distribution will be peaked 
around photon numbers where the product of the cosinusoidal and sinusoidal distributions 
is maximised. Examination of figures 1(a) and 2(a) shows that this will be around five and 
ten photons. This is borne out by figure 4 which shows two typical results for six detected 
atoms, alternating between detection in the excited and ground states. The only significant 
probabilities are for 4-5 and 9-12 photons. Other orderings give slightly different ratios of 
probability between these numbers, but the results are similar. 

V. CONCLUSIONS 

In this paper we have analysed field measurements in the micromaser. The field is not 
directly measurable, but knowledge of the prepared initial atomic state and the measured 
state after the interaction between the atoms and the cavity gives information about the 
initial cavity photon number distribution. If nothing is known a priori about the cavity 
field this allows us to calculate the POM elements for the field corresponding to particular 
sequences of atomic detections in either the excited or ground states. 

We have considered the case where the atoms are all initially prepared in the excited 



13 



state. If an atom is detected in the excited state the calculated initial photon number 
distribution is simply the distribution after the measurement multiplied by a cosinusoidal 
Rabi factor which depends upon one more than the final photon number. This happens for 
each repeated detection in the excited state. Atoms detected in the ground state give the 
initial field a sinusoidal distribution which depends upon the one more than the final photon 
number. 

If the atoms are measured in the excited state the cavity is much more likely to have 
been in a photon number state corresponding to a trapping state, for which the Rabi phase 
shift is a multiple of ir. Conversely, if the atom is measured in the ground state the cavity 
field cannot have been in a trapping state. The probability of states with these photon 
numbers is zero. Repeated measurement in the ground state therefore measures photon 
number distributions with large gaps. This can be used both to measure and prepare photon 
number states, particularly where some information about the field also exists. 

Earlier work on retrodictive master equations 0, and on the derivation of the POM 
elements corresponding to homodyne detection (the line states) [T3|], has shown that the 
boundaries between preparation, evolution and measurement in quantum theory are arbi- 
trary choices. This is emphasised in the micromaser where the measurement and preparation 
of particular atomic states corresponds to a single measurement of the field. A measurement 
of a quantum system necessarily includes all of the available information about the system. 
If any of this information should change so does the measurement POM. If the atoms had 
been prepared in a superposition of the excited and ground states then measurement in the 
excited or ground state would correspond to different field POM elements to those described 
here. We will explore preparation and measurement of the atom in superpositions in further 
work. 
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APPENDIX 

The sequence of detections of atoms in their ground or excited states represents a mea- 
surement of the initial field. Our analysis has concentrated on the retrodictive field state 
prior to interaction with the first atom. It is not necessary, however, to introduce the retro- 
dictive state in order to analyse the field measurement. It can be calculated simply by 
normalising the derived POM element. In calculating the POM elements we have placed 
a limit on the photon number of the final cavity field. This is also not strictly necessary 
provided that we are only interested in the field measurement and its associated POM. If 
we know the initial field state, Pj Te< ^(tp), then these POM elements give us the probabilities 
for any given sequence of atom detections. 

In constructing the POM, we start with the unmeasured field state after the final (s th ) 
atom has passed throught the cavity. We associate with this the trivial POM element 

oo 

flf^ = l { =J2\n)(n\. (A.l) 

n=0 

The POM element associated with immediately prior to detection of the atom (denoted 
n|) will depend upon whether the atom was detected in the excited or ground state: 

Bf(e) = (e|f> t |e)(e| <8> l { U\e) 

= cos [n(n + l)r/2] l f cos [Q(n + l)r/2] 

= cos 2 [Q(n + l)r/2] , (A.2) 
15 



fq(g) = (e\Ui\g)(g\®i i U\e) 

= sin [Q(n + l)r/2] El { E ] sin [fi(n + l)r/2] 

= sin 2 [ft(n + l)r/2]. (A.3) 

Here n is the photon number operator and E(E^) is the bare lowering (raising) operator 







E = J2\n){n+l\ (A.4) 

n=0 

oo 

£t = £|n+l)(n|. (A.5) 



n=0 

These operators have the properties 

Ef(n) = f(n + 1) (A.6) 

/(n) = #/(n + l) (A.7) 

EE^ = 1 (A.8) 

£^£ = 1- |0> <0|. (A.9) 

We can use these operators to write the POM elements associated with any given se- 
quence of atom detections by induction. This means starting with the last atom detected 
and writing the POM element cos 2 [0(n + l)r/2] if the atom was in its excited state and 
sin 2 [Q(h + 1)t/2] if it was in its ground state. For each preceding detection we pre- and 
postmultiply the POM element by cos [Q(n + 1)t/2] if the atom was measured in its excited 
state. If the atom was found in its ground state then we premultiply by sin [Q(n + 1)t/2] E 
and postmultiply by E^ sin [Q(h + l)r/2]. 

As an example, we consider a sequence of three atoms. The corresponding POM elements 
lif(z, j, k) (i, j, k = e,g) are associated with the atoms interacting with the cavity in the order 
i,j then k. The eight possible POM elements are. 

fl f (e,e,e) = c 6 (l), (A.10) 
n f (e, e, g) = c 2 (l)S(l)££tS(l)c 2 (l) = c 4 (l)S 2 (l), (A.ll) 
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TT / \ 

nf(e,^,e) = 


c(l)s(l)Ec (l)E T s(l)c(l) = c (l)s (l)c (2), 


/A -1 r\\ 

(A.12) 




c(l)s(l)Es(l)EE T s(l)i? T s(l)c(l) = c (l)s (l)s (2), 


(A.13) 


TT / \ 

Hfl^, e, e) = 


s(l)Ec (l)_E T s(l) = s (l)c (2), 


(A.14) 


rif(#,e,#) = 


5(l)£c(l)5(l)£#5(l)c(l)£ t 5(l) = s 2 (l)c 2 (2)s 2 (2), 


(A.15) 


n f (#,#,e) = 


s(l) J E5(l) J Bc 2 (l) J B t s(l) J Ets(l) = s 2 (l)s 2 (2)c 2 (3), 


(A.16) 


tif(g,g,g) = 


5(l)£5(l)£5(l)££ t 5(l)£ t 5(l)£ t 5(l) = s 2 (l)s 2 (2)s 2 (3). 


(A.17) 



Here we have written 

c(m) = cos [Q(h + m)r/2] , (A. 18) 

s(m) = sin [ft(n + m)r/2] . (A.19) 

We can also summarise the construction of the field POM elements by iteration starting 
from the first atom to enter the cavity: 

(1) Start with the first atom to interact with the cavity. 

(2) If the atom is measured in the excited state then write c 2 (l). If the atom is measured 
in the ground state them write s 2 (l) and add 1 to the arguments of all subsequent c and 5 
operators. 

(3) Follow the same prescription for each subsequent atom in turn. 
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Figure Captions 

Figure 1. Measured cavity photon number probability distributions for atoms detected in 
the excited state (a) one atom, (b) five atoms. 

Figure 2. Measured cavity photon number distributions for atoms detected in the ground 
state (a) one atom, (b) three atoms, (c) six atoms. 

Figure 3. Measured cavity photon number probability distributions for two atoms detected 
in the ground state given that there were no more than three photons in the cavity prior to 
the passage of the atoms. 

Figure 4. Measured cavity photon number probability distributions for six detected atoms, 
three in the ground state, three in the excited state in the order (a) g, e, g, e, g, e and (b) 
e, g, e, g, e, g. 
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